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\ This is the second part of the series of papers on symmetry properties of a class of vari- 

■ able coefficient (l+l)-dimensional nonlinear diffusion-convection equations of general form 
f{x)ut = {g{x)A{u)ux)x + h{x)B(u)ux- At first, we review the results of [12] on equivalence 
transformations and group classification of the class under consideration. Investigation of 

f~| I non-trivial limits of parameterized subclasses of equations from the given class, which gcner- 

C^H' ate contractions of the corresponding maximal Lie invariance algebras, leads to the natural 

fH \ notion of contractions of systems of differential equations. After a brief discussion on con- 

tractions of symmetries, equations and solutions in general case, such types of contractions 
are studied for diffusion-convection equations. A detailed symmetry analysis of an inter- 
esting equation from the class under consideration is performed. Exact solutions of some 
subclasses of the considered class are also given. 
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The present manuscript continues the series of works on symmetry properties of the nonlinear 
variable-coefficient diffusion-convection equations of form 

f{x)ut = {g{x)A{u)ux)x + h{x)B{u)ux, (1) 
where / = /(x), g = g{x), h = h{x), A = A{u) and B = B{u) are arbitrary smooth functions of 



^ ■ their variables, f{x)g{x)A{u)^^. 



Based on the results on equivalence transformations and group classification adduced in 
the first part [12] of the presented work, we continue studying Lie group properties of the 
equations under consideration. Namely, investigation of non-trivial limits of parameterized 
subclasses of equations from class (1), which generate contractions of the corresponding maximal 
Lie invariance algebras, leads us to the natural notion of contraction of (systems of) differential 
equations. 

It is well-known that exponential cases of parameter-functions, which admit extensions of 
maximal Lie invariance algebras, are often limits of the power ones with arbitrary exponents [1, 
2,25]. By analogy with terminology accepted for the Lie algebras we will call such limits as 
contractions. In such situations some authors exclude exponential cases from formulation of 
final results of group classification that is correct only under explicit admission of combined 
usage of both point equivalence transformations and contractions in the framework of group 
classification. Another famous contraction is one of the 1-dimensional Ricci flow ut = A Inn 
from the porous medium equation ut = m~^A(n"^ — 1) under the limit m ^ [28]. 
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All the above-mentioned contractions were found by ad hoc procedures. To the best of our 
knowledge this is the first work giving the precise definition and mathematical background for 
contractions of equations, algebras of symmetries and solutions. The value of the presented 
theory is illustrated in Section 3 using as example the variable-coefficient diffusion-convection 
equations (1) investigated. 

Such contractions allow us to establish additional connections between different cases of 
extension of maximal Lie invariance algebra in class (1) and ones between solutions of equations 
in these cases. Thus, e.g., applying contractions to the known solutions of initial equations one 
can easily obtain exact solutions of the target equations. 

In fact, group analysis (in particular, invariant and partially invariant solutions and solutions 
obtained from separation of variables) is the only systematic method we know for deducing exact 
solutions of nonlinear partial differential equations, and the role of exact solutions in our un- 
derstanding the mechanics and mathematics of mathematical models cannot be over-estimated. 
Classes of invariant solutions usually include self-similar (automodel) solutions arising from scal- 
ing invariance, travelling waves which are invariant with respect to translation symmetries and 
other classes of solutions having physical and analytical importance. For most of the nonlinear 
systems invariant solutions are the only available exact solutions. 

Below we construct exact solutions of equations from class (1) using different approaches. 
We start by illustration of the simplest approach for deriving exact solutions of the complicated 
equations from the ones of a simpler model by means of application of equivalence transforma- 
tions. At first sight such investigation may seem to be not very interesting and trivial. However, 
a number of authors investigate symmetry properties of different classes of equivalent equations. 
In such way they perform a huge number of unnecessary cumbersome calculations. Since equiv- 
alence transformations may be quite complicated, sometimes it seems to be impossible to obtain 
directly complete and correct results for some cases which are equivalent to simple ones, although 
these results can be easily reconstructed with application of equivalence transformations. 

At the same time, it is a general mathematical rule that equivalent in some sense objects pos- 
sess equivalent properties. In particular, if two systems of equations are equivalent with respect 
to point transformations then there exists a one-to-one correspondence between their maximal 
Lie invariance algebras, spaces of conservation laws, potential symmetries, exact solutions etc. 
All the above mentioned features for more complicated models can be constructed from ones 
of the simpler equivalent models by means of application of known equivalence transformation 
easier than by direct calculations. 

We performed the classical Lie reduction procedure for an equation that cannot be reduced 
to a constant coefficient equation from class (1). Then we investigate in more detail an equation 
having, from the one side, remarkable Lie symmetry property (namely, s/(2, M)-invariancc) and, 
from the other side, admitting application of powerful non-Lie methods of construction of exact 
solutions. 

The present paper is organized as follows. First of all, in Section 2, for convenience of the 
reader we review (in the most suitable for our purpose form) the results of [12] on equivalence 
transformations and group classification of equations (1) using simultaneously two different 
gauges g = 1 and g = h. Then, in Section 3 we introduce the notion of contractions of equa- 
tions and contractions of symmetry algebras, give a necessary theoretical background and found 
contractions of previously classified diffusion-convection equations and those of corresponding 
symmetry algebras. Examples of application of additional equivalence transformations are shown 
in Section 4, where we reconstruct exact solutions of variable coefficient equations with four- 
dimensional symmetry algebras (reducible to constant coefficient ones) from the known solutions 
of constant coefficient equations. A case of equation reducible to a constant coefficient form that 
do not belong to the class under consideration is presented in Section 5. In Section 6 an en- 
hanced group analysis of essentially variable coefficient sZ(2, M)-invariant diffusion-convection 
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equation is presented and a set of exact Lie and non-Lie solution is constructed via reduction 
with respect to the classical Lie and nonclassical symmetries. 

Results of Section 3 will be generalized in the third part [13] of this series of papers where we 
introduce the notion of contractions of conservation laws, study the space of local conservation 
laws and investigate contractions of the conservation laws of equations from class (1). Ibid 
wc introduce also some generalizations of equivalences of conservation laws that allow us to 
generalize procedure of construction of potential systems for systems of differential equations. 
Using this new technique we construct all possible inequivalent potential systems for equations 
from class (1). 

2 Equivalence transformations and group classification 
of diffusion— convection equations 

In this section we briefly summarize the results of the first part [12] of this series of papers 
concerning equivalence transformations and group classification of class (1), presenting these 
results in the most suitable form for our purpose. Due to lack of power of the usual equivalence 
group G~ it seems to be impossible to obtain result of group classification of class (1) up to 
in explicit form [15]. So, we adduce the extended equivalence group [11,12] containing transfor- 
mations depending nonlocally on arbitrary elements. Using the direct method we construct the 
complete in this sense group G~ formed by the transformations 

t = 5it + 62, x = X{x), u = Ssu + S4, 

f = A^f, g = eie^-^X^ipg, h = £ie^^ifh, A = e2A, B = £s{B + e4,A), 

where 6j [j = 1,4) and Si {i = 1,4) are arbitrary constants, diSs£i£2S3 ^ 0, X is an arbitrary 

smooth function oi x, X^ ^ 0, (p = e ■' . 

Group G~ contains a normal subgroup of gauge equivalence transformations [11,12,18] 

f = £iipf, g = £i£-^ipg, h = £i£^'^Lph, i = £2-4, B = £^{B + £4^A), (2) 

where (p = e ai^) ^ {j, = are arbitrary constants, £i£2S'i 7^ (the variables t, x and 
u do not transform!). The transformations (2) act only on arbitrary elements and do not really 
change equations. Application of gauge equivalence transformations is equivalent to rewriting 
equations in another form. In spite of really equivalence transformations, their role in group 
classification comes not to choice of representatives in equivalence classes but to choice of form 
of these representatives. 

Choice of a gauge for the arbitrary elements is very important for solving and for the final 
presentation of results. It is more convenient to constrain the parameter-function g instead of / 
in class (1). The next problem is choice between gauges of g. Group classification of (1) is 
performed in [12] in two different gauges: g = 1 and g = h. For further investigation wc should 
reduce these results to the floating gauge in order to obtain the simplest representatives in each 
equivalence class. So, below, in contrast to [12] we adduce the group classification with floating 
gauge. The cases B<^{l,A) and A = I are investigated in the gauge g = h. In the other cases, 
to obtain results in a simpler explicit form one should use the gauge g = 1- 

Theorem 1. The Lie algebra of the kernel of principal groups of (1) is A'^ = (dt). A complete 
set of G'^ -inequivalent equations (1) which have the wider Lie invariance algebras than A'^ is 
exhausted by cases given in tables 1-3. 
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In tables 1-3 we list all possible G~-inequivalent sets of functions f{x), h{x) A{u), B{u) 
and corresponding invariance algebras under the gauge g = 1 or g = h. Numbers with the 
same Arabic numerals and different Roman letters correspond to cases that are equivalent with 
respect to additional equivalence transformations. Explicit formulas for these transformations 
can be found in [12]. The cases which are contained in tables with different Arabic numbers or 
are numbered with different Arabic numerals are reciprocally inequivalent with respect to point 
transformations, except the cases with asterisked numeration, presented here to preserve the 
same numeration as in [12]. 

The operators from tables 1-3 form bases of the maximal invariance algebras iff the corre- 
sponding sets of the functions /, h, A and B are G'^^'-inequivalent to ones with more extensive 
invariance algebras. For example, in case 3.1 the adduced operators have the above property iff 

Case 1.2a' is equivalent to case 1.2a with respect to transformation 
i=t, x = ln|x|, u = u, A = A, B = B — A, p = p + 2 
from We adduce case 1.2a' here for the convenience of presentation of results only. 

Table 1. Case of V^(m) (gauge g = 1) 



N 


B{u) 


fix) 


h{x) 


Basis of A'"^" 


1 


V 


V 


V 


dt 


2a 


V 




1 


dt, ptdt + dx 


2a' 


V 




X-' 


dt, {p + 2)tdt + xdx 


2b 


1 






dt, e-\dt-dx) 


2c 


1 




x\x\P + (ix-^ 


dt, e-'-p+''^\dt - xdx) 


3 


1 


-2 

X 


x^^ In \x\ 


dt, e~*xdx 


4a 





1 


1 


dt, dx, 2tdt + xdx 


4b 


1 


1 


x 


dt, e-'dx, e-^\dt-xdx) 



Here p 6 {0, 1} mod Gi in case 2; p ^ —2 in case 2c; /3 € {0, ±1} in case 2b. 



Table 2. Case of A{u) = e*"" 



N 


B{u) 


fix) 


gix) 


h{x) 


Basis of A™*'' 


1 





V 


1 


1 


dt, tdt - du 


2 


e"" 




1 


\xV 


dt, {p^l - pu — 2u - qfi + fi)tdt + ip — v)xdx + (? + 


2* 


e"" 




1 




dt, (pM -pv - n)tdt + (m - '^)dx + du 


3 


we" 






gP^" 


dt, {2p + q)tdt + dx - 2pdu 


4 




1 


1 


1 


dt, dx, ip — 2u)tdt + Hp - v)x + i'xt)dx + du 


5 


u 


1 


1 


1 


dt, dx, tdt + ix- t)dx + du 


6a 





fix) 


1 


1 


dt, tdt - du, atdt + iPx^ + 71a; + io)dx + fixdu 


6b 


1 


\x\^ 


1 


ex\x\P 


dt, xdx + (p + 2)du, e-^(P+2)*(at - exdx) 


6b* 


1 




1 


ee" 


dt, dx + du, e~^\dt - edx) 


6c 


1 


x-^ 


1 


ex^^ 


dt , xdx , tdt — stxdx — du 


7a 





1 


1 


1 


dt, tdt - du, 2tdt + xdx, dx 


7b 


1 


1 


1 


1 


dt, dx, tdt — tdx — du, 2tdt + (x - t)dx 


7c 


1 


1 


1 


ex 


dt, xdx + 2du, e-'*dx, e-^^\dt - exdx) 


7d 





X-' 


1 


1 


dt, tdt - du, xdx - du, x^dx + xdu 


7c 


1 




1 


x-^ 


dt, xdx - du. e'idt - xdx), e\x^dx + xdu) 



Here (/x, f) £ {(0, 1), (1, i>)}, v ^ jj, in cases 2, 2* and 4; = 1 and 1 in the other cases; g ^ — 1 in 

case 2* (otherwise it is subcase of the case 1.2a'); e = ±1 in cases 2, 6b-6c and 7e; p ^ {—3, —2,0} in case 6b; 
a, /3, 71, 70 = const and 

-Wx - 271 + Q 



/ [x] = exp 



Px'^ + 71a; + 70 
Case 2{q = —1) is a subcase of case 1.2a' 



dx 
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Table 3. Case of A{u) = |u|' 











9^x) 




IDdoiD Ui r\. 


1 


V 





V 


1 


1 


dt, iJ.tdt - udu 


o 

Z 


V 


Pl 


Fl 


1 
1 


\X\ 


dt, in + PH — — pi' — 2v)tdt 
— i/)a;dj: + (g + Ijudu 


9* 


w 

V 


|M| 




1 
± 


c 




Q 
O 


V 






C 


„P^^ 


Ot, {^fJ'P "T Q)iOt Ox zpuOu 


4 


V 


1 




1 


exf{x) 


dt. 





n 
U 


w 

V 


e 


(-. 


e 


Ot, e Ox 


D 


n 
u 


w 

V 


e 


c 


e 


Ot, 6 yOt — 7Cx } 


7 


n 


It 


t. 




c 


Ct, lut ^P^u 


8 


V 




1 


1 


1 


dt, 

{ II 9;/^fr)^ -\- ( ( II 7/W -1- t/-i^i\fi -X. iifi 


Q 


w 

V 


ill tl 


J. 


1 

L 


1 

J. 


Ot, <Jx, t^iut yH'X — iJOx uuu 


iU 


u 


It 


6 


C 


e 


Ot, e Ox, Ox — ^pOu 


il 


U 


In L.I 
111 \U\ 






«p^c^ 


yt, e Ox, Ox — 2puOu 


12a 


V 







1 


1 


dt, fj.tdt - itc?u, 

aiot + [(n + ^)px + 71a: + 'YojOx + pxuOu 


1 9K 
iZD 


w 
V 


1 
1 


|x| 


1 
1 


IP 


fl ii'yfl -i- (n J- 0\',if) o~^(P+^)*ffl c-ffi \ 

Ot, uxOx + \JP + 'ijuOu, e (Ot — sxOxj 


Izu 




1 

i 


a: 

e 


1 
i 


ee 


fl ,,fi 1 r,,fl ^,^^^tfl ^fl ^ 

Ot , i-iOx + uOu , e {Ot — sOx ) 


IZC 




1 

± 




1 

± 






13 


—6/5 


1 




1 




t^Q^ + ^2tx + x^)dx — 5{t + x)udu 


14a 


/-4/3 





1 


1 


1 


dt, fJ,tdt - udu, dx, 2tdt + xdx 


14b 


/-4/3 


1 


1 


1 


1 


dt, )J,tdt - ^J,tdx - udu, dx, 2tdt + {x — t)dx 


14c 


7^-4/3 


1 


1 


1 


ex 


dt, nxdx + 2udu, e-^'dx, e-''^\dt - exdx) 


14d 


^ -4/3, -1 





|a;| ''+1 


1 


1 


dt, i^tdt - udu, in + 2)tdt - (/i + Vjxdx, 








|a;| f'+i 




ea;|a;| '■+1 


(/i + l)x^dx + xudu 


14e 


-4/3, -1 


1 


1 


dt, m(m + '^)xdx — + 2)udu, 

e ''+1 \at—exOx),e ((/Lt+ijx Ox-yxuOu) 


14f 


-1 







1 


1 


dt, tdt + udu, dx — udu, 2tdt + xdx — xudu 


14g 


-1 

— i 


-1 

1 


e 


-1 
i 


£6 


Ot, Ox — UOu, (X + et — 2)Ox — [X + etjuOu, 

e-''{dt-edx) 


1 /I Vi 


— Z 


1 
i 




1 
i 




f) ryfl + fl Oc + ry'fl _L ^ i fl 

Ot, xOx, ztot — zstxOx + UOu, 
e^*{x^dx — xudu) 


ioa 


-4/3 


U 


-1 
1 


1 
i 


1 


Ot, 'aot + ouOu, Ox, Mot + XOx, 
x^dx — .ixudu 


iOD 




1 
1 


1 
1 


1 
1 


1 
1 


Ot, ^lOt -r ^XOx — oUOu, ^lOt -]- \X — tjOx, 

dx, {x + t^dx — 3(x + t)uBu 


15c 


-4/3 


1 


1 


1 


ex 


dt, 2xdx — Sudu, e~^*dx, 

e-^^\dt - exdx), e^'ix'^dx - Zxudu) 


16 





u 


1 


1 


1 


dt, dx, tdx - du, 2tdt + xdx - udu, 
t^dt + txdx — {tu + x)du 



Here v ^ \x; e = ±1; g ^ — 1 in case 2* (otherwise it is subcase of the case 1.2a'); p 7^ —2, — (S/i + 4)/(p + 1) in 
case 12c; a, /3, 71, 70 = const, and 



/ (x) = exp 



-(3m + 4)/3x-3 
+ l)^a;2 + x 



dx 



fix) 



exp 



-(3/i + 4)/3x-27i+Q 
{fl + l)(3x^ + 71a; + 70 



dx 
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3 Contractions of equations and of Lie invariance algebras 



In this section wc investigate non-trivial limits of parameterized subclasses of equations from 
class (1), which generate contractions of the corresponding maximal Lie invariance algebras 
as realizations of abstract Lie algebras in the complete space of (independent and dependent) 
variables. By analogy with terminology accepted for abstract Lie algebras, we will call such 
limits of parameterized classes and their Lie symmetries as contractions of equations and Lie 
invariance algebras correspondingly. 

First we introduce some related notions and formulate two simple but very useful statements 
concerning contractions of equations in the general case. 

Consider the class {C^} of systems C^: L{x, U(^p^, A) = of ^ differential equations for m un- 
known functions u = {u^, . . . , u^) of n independent variables x = {xi, . . . , which are param- 
eterized with the parameter A. Here denotes the set of all derivatives of u with respect to 
X of order not greater than p, including u as the derivatives of the zero order. L = {L^, . . . , L') 
is a tuple of I fixed functions depending on x, U(^p^ and A. For simplicity we assume A as a 
single numeric (real or complex) parameter. (Extension to more general case with respect to A 
is obvious.) We also suppose that the systems are totally nondcgcncrate. 

Let jC'^ = {L^{x,U(^p^, X) =0, k = 1, ... ,1} he a maximal set of algebraically independent 
differential consequences of that have, as differential equations, orders not greater than p. 
All such sets of differential consequences determine the same manifold in the jet space 
Let us fix a value of the parameter A = Aq and a point Jq G jC^° . 

Definition 1. The systems jC^ weakly converge to the system jC^° near the point jo € J0-^° 
under A — >^ Aq, iff for any A from a deleted neighborhood A of Aq there exist and an open 
neighborhood of jo in J^^^ such that 

1) L^{x, U(p-j , X). k = 1, . . . ,1 and their partial derivatives with respect to x and 7i(pj converge 
to differential consequences of and their corresponding derivatives pointwise on under 

2) if L^, k = 1, . . . ,1 are extended to Aq by continuity then for an /-element subset of the 
variables the Jacobian \dL^/dv^ \ ^ on Q for any A from a neighborhood of Ao. 

We will use the denotation C'^", A ^ Ao. 

Consider the parameterized set {Q^, A G A} of the operators 

n m 

i=l j=l 

on the space of variables {x, u) and a fixed operator 

n m 

Q' = Y1 ^°'(^' + E ^)^«^- • 
i=i j=i 

The notation pr^^) Q^: ^ ^ denotes that the coefficients of the standard j3-order 

prolongation of to J^^ pointwise converge to the corresponding coefficients of the p-order 
prolongation of the operator on the set $7. 

Lemma 1. Let for any X from a deleted neighborhood of Xq the system be Lie invariant with 
respect to an operator Q^, — ^ and pr^p) — ^ P^(p) ' A — ^ Ao- Then the system is 
Lie invariant with respect to the operator . 

Proof directly follows from the infinitesimal invariance criterion. 
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Corollary 1. Suppose that for any X from a deleted neighborhood of Xq the system is Lie 
invariant with respect to a fixed operator and with respect to an operator sufficiently 
smoothly depending on value of X, andQ'^ are linearly independent, C'^ C'^" andpr^^p^ Q'^ 
pr^p) Q*^, A ^ Aq. Then the system C^^ is Lie invariant with respect to the operators and 
d'^Q^/dA^I^^^^. Here d^Q^/dX^ is the operator with coefficients being k-th order derivatives of 

the corresponding coefficients of the operator with respect to X, and k is the first integer for 
which d''Q^/dX''\^^^^ and Q'^ are linearly independent. 

Proof. Both pairs = = and 



A=Ao 



satisfy the conditions of the previous lemma. □ 

Note 1. For contracted systems of differential equations we can introduce the notion of con- 
tractions of solutions and then investigate connections between contractions of Lie invariant 
solutions and contractions of corresponding subalgebras of the Lie invariance algebras (sec Sec- 
tion 8). It is also possible to consider contractions of involutive sets of reduction operators and 
solutions which are invariant in the non-classical (conditional) sense. The study of the above 
contractions and contractions of symmetries and invariant solutions of other kinds will be the 
subject of our further investigations. 

There exists a number of different non-trivial contractions of equations from class (1) and, 
in particular, equations having wider Lie invariance algebras than the kernel algebra A'^ = (dt). 
Contractions allow us to find additional connections between different cases of extension of 
maximal Lie invariance algebra in class (1) and between solutions of equations in these cases. 

The simplest contractions are the ones in parameterized subclasses of equations with respect 
to their parameters. For instance, case 2.6c is the limit of case 2.6b under p ^ —2. The Lie 
invariance algebra of case 2.6c can be obtained from that of case 2.6b by means of the contraction 
with p — ^ —2: 

e{p + 2) 

Hereafter the operators are taken from the corresponding tables. The "left-hand-side" ("right- 
hand-side") operators are from the given basis of the case which is the object (target) of the 
contraction. The basic operators are numerated in accordance with their order at the tables. 

The above contraction illustrates corollary 1 since Q^^ + eQ^^ Q^^ = Q^, p ^ —2 and 
Qi = d(Q^P + eQ'^P)/dp\ The tuples (Q^p , Q'^p , Q^p) and ((3^Q^Q^) determine equivalent 

third-rank realizations of the algebra A2.1 © Ai. We use the Mubarakzyanov's notations [19] of 
low-dimensional solvable real Lie algebras. Ai denotes the one-dimensional Lie algebra and A2,i 
does the non-Abelian two-dimensional Lie algebra. It follows from results of [23] that up to weak 
equivalence of realizations there exists the unique third-rank realization of the algebra ^42.1 ©^1. 
Therefore, equivalence of realizations {Q^p , Q^p , Q^p) and {Q^,Q^,Q'^) is foreknown. 

In table 3 an analogous contraction is 12.b — ^ 12. c, p —2, where 

£{p + 2) 

Existence of analogous cases of extensions of maximal Lie invariance algebras and analogous 
contractions in tables 2 and 3 can be explained in terms of contractions. 
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Consider equations of the form 

f{x)ut = ig{x)\u\''ux)x + h{x)B{u)ux, 

covering the cases from table 3 (or 3'), which are parameterized by /i. The parameter-function B 
takes the values 

Inl*^, \u\'^ + X, 1, 0, I'ul'^lnlnl, Inl-ul. 

We do not specify here the values of f{x) and h{x) explicitly. Limits with respect to exponents 
of \u\ are more cumbersome than the already considered ones. At first we have to carry out the 
transformation parameterized with a parameter 5: 

u 

u = l + -, ii = 5fl, i' = 5v, a = 5a, 71 = (57i, 70 = ^7o- 


Hereafter we change parameters iff they are in the corresponding equation. For two latter values 
of B we have additionally to transform the independent variables and the parameters p and q: 

t = S'^i, X = Sx, p = S'^p, q = Sq. 

Then, we proceed to the limit S — > +00. The limit values are A = e''" and 

e^", e^" + x, 1, 0, iie'^", u 

for B (correspondingly to the order of the above values of B). For the most of contracted cases 
the parameter- functions / and h do not formally change their values. The exceptions are the 
cases 3.12a, 3.14d and 3.14e where / and h depend explicitly on fi. The described way results 
in the following contractions: 

3.1^ ^ 2.1, 3.2,^1/pq > 2.2^j>pg, 3.2 ^jyp > 2.2 jicpj 3.3^pq > 2.3pg, 3.8^;/^^ > 2.4j>pq, 
3.9^^2.5, 3.12a^„^^j^„ ^ 2.6aQ,/3^^^o, 3.12bp ^ 2.6bp, 3.12b* ^ 2.6b*, 
3.14a^ 2.7a, 3.14b^ 2.7b, 3.14c/, ^ 2.7c, 3.14d^ 2.7d, 3.14e^ 2.7e. 

Here /i = 1 in all the above cases from table 2 excluding case 2.2. 

Analogous contractions can be performed in cases when the parameter-functions / and h are 
powers by means of the transformation 

X 1 1 
t = (5^~^^t, X = 1 -|- -r, p = Sp, q = Sq, u = u -\ In S (or u = S"-^u) 

u — fj, ^ ' 

for cases from tables 2 (or 3 correspondingly) and proceeding to the limit b — > -|-oo: 

^•^iwpq ^ 2.2 ^i/p,q=ij 2.6bp > 2.6b ^=1, 3.2^^^^ > 3.2 ^up,q=\i 3.12bp > 3.12b . 

To obtain similar contractions of cases from table 1, we have to use another transformation of 
the variables and parameters: 

i X r A ^ 

t=-z, x = l + -, u = u, p = op, A = — 
00 

that leads to the contractions 1.2a'pA 1.2a-^, 1.2cp^^ — > 1.2b-^^. 

The question whether the above contractions exhaust non-trivial ones between equations 
from class (1), which admit extensions of maximal invariance algebras, remains open. 
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4 Construction of exact solutions 
via equivalence transformations 

Results of group classification and additional equivalence transformations can be used for con- 
struction of exact solutions of equations from class (1) admitting extensions of maximal Lie 
invariance group. In [15,24] known exact solutions of "constant coefficient" diffusion-convection 
equations were mapped to ones of variable coefficient equations (1) by means of additional 
equivalence transformations. 

In view of Theorem 5 of [12] if an equation of form (1) is invariant with respect to a Lie 
algebra of dimension not less than 4 then it can be reduced by point transformations to one 
with f = g = h = 1. Therefore, we can reconstruct solutions of these variable-coefficient 
equations from the already known solutions using additional equivalence transformations found 
in [12]. Thus, e.g., equation 3.14g 

e^«t = {u~^Ux)x + ee^u^ (3) 

can be transformed to the fast diffusion equation ut = {u~^Ux)x with the point transformation 

i = e^*/e, x = x + et, u = e^+^*u. (4) 

The list of known G^'^^-inequivalent Lie solutions of the fast diffusion equation is exhausted 
by the following ones [26]: 



2t -2t 2t 2sin2i 

5j n = b) u = 7) u = — — °} 



cos^ X 



cosh^ X ' sinh^ x ' cos 2t - cos 2x ' 



2sinh2t , 2sinh2t 

9)u = — — — , 10) u = -- 



(5) 



cosh 2x — cosh 2t ' cosh 2x + cosh 2t ' 

2cosh2t , 2sin2t , 2sinh2t 

11)^=^171^ ^ITTT' 12) u = 7— —, 13) u 



sinh 2x — sinh 2t ' cosh 2x — cos 2t ' cosh 2t — cos 2x ' 

Here S = ±1, C is an arbitrary constant; 1)— 7) are Lie invariant solutions and 8)— 13) are obtained 
with nonclassical potential symmetries. 

Applying transformation (4) to (5) we obtain a list of exact solutions of the variable-coefficient 
equation (3): 



e 



-{x+£t) g.g-(a;+et) 



u = - : — u = 



1 -I- ,5e^+£*+e^*/£ ' sx + eH - e^* -I- C'e^*-^(^+^*)^" 

2c-''' 2c-'' 

U = T, U 



£{x + et)"^ + 6e^^^' £cos^{x + sty e cosh^ (x -I- et) ' 

2e-^ _ 2e-(^+^*) sin(2e^Ve) 



u = 



£:sin^(x -|- et) cos(2e^*/^) ~ cos 2{x + et) 

2e-(^+^*) sinh(2e^V£) ^_ 2e-(=^+^*) sinh(2e^7e) 



cosh2(x -I- et) - cosh(2e^V£) ' cosh2(a; -|- et) + cosh(2e^V£) ' 

2e-(^+^*) cosh(2e^7e) ^_ 2e-(^+^*) sin(2e^V£) 



u = 



sinh2(a; -|- et) - sinh(2e^V£) ' cosh2(a; -|- et) - cos(2e^V£) ' 

2g-(x+£t) sinh(2e"Ve) 
cosh(2e^*/£) — cos 2{x + et) ' 
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The lists of known exact solutions and Lie reductions of the constant coefficient and some of 
the variable coefficient equations can be found, e.g., in [10,21]. Using them one can easily obtain 
exact solutions of any of the variable coefficient equations having the Lie symmetry algebra of 
dimension not less then 4. Since the solutions of these equations can be reconstructed from ones 
presented in other sections, we can turn back to the more interesting cases. 

5 Example of Lie reduction 

In this section we consider in more details Lie reductions of variable coefficient equation 3.10 

ef-^t = (ef"\,), + ef^'^xtx„ (6) 
which is invariant with respect to the three-dimensional Lie symmetry algebra 

{dt, e-^P'd,,d,-2pdu). 

In contrast to the case of equations with four-dimensional Lie symmetry algebra wc cannot 
reduce equation (6) to an equation of form (1) with f = g = h = 1. However, it is an interesting 
feature of this equation that using a point transformation v = u + 2px it can be mapped to a 
constant coefficient reaction-convection-diffusion equation 

vt = Vxx + vvx - 2pv (7) 

that does not belong to class (1). Similarly to Section 4 for simplification of the technical 
calculations we will investigate the constant coefficient equation (7) instead of (6). The Lie 
symmetry algebra of equation (7) 

(Xi = dt, X2 = e-^P\d^ + 2pd^), Xs = d^) 

is a realization of A2.1 © Ai [19]. These operators generate the following group of point trans- 
formations: 

i=t + ei, 5 = x + £2e'~^^* + £3, V = V + 2e2pe''^P^ . 
A list of proper inequivalent subalgebras of the given algebra is exhausted by the following ones 

where a and are arbitrary constants, e = 0, ±1 [20]. 

The first three (one-dimensional) subalgebras lead to Lie reductions to ordinary differential 
equations, the fourth and sixth (two-dimensional) ones yield reductions to algebraic equations. 
Lie reductions with respect to these subalgebras are summarized in table 4. One can easily 
check that it is impossible to construct a Lie ansatz corresponding to the subalgebra (X2, X^). 



Table 4. Lie reductions of equation (7). 



N 


Subalgebra 


Ansatz V = 




Reduced equation 


1 
2 
3 
4 
5 


(Xi + aXs) 
(X3 + eX2) 

{X2) 
{Xi, Xs) 
(Xi-F/3X3, X2) 


tp{(jj) + 2px 
C 

2px - 2pPt + C 


X — at 
t 
t 


(p" + {v> + a)(p' - 2pip = 
2pe'^'^ 

(p' = 

c = o 

-2pP = 
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We are succeeded in solving equations 2-5 from table 4 that give us the following invariant 
solutions of equation (7): 



V = 0, V = 2px + C, V 



2pe.r + Cf:2p* 



g2pt _|_ g. 



The corresponding exact invariant solutions of equation (6) have the form 



u = —2px, u = C, u = 



2psx + Ce^P* 



(8) 



g2pt _|_ £ 



where C is an arbitrary constant. 

Ansatzes 4.2 and 4.3 give a hint for a possible form 

V = ip{t)x + V'(i) 

of nonlinear separation of variables for construction of exact solutions of equation (7). Substi- 
tution of the ansatz to equation (7) leads to antireduction: 

(p' = (p'^ — 2pip, tp' = (ftp — 2p^. 

Solving the above system of ODEs for (p and ip we obtain exactly the solutions of equations 4.2 



Note 2. Using the point transformation i = e x = x,v = e^^^v equation (7) can be mapped 
to a variable coefficient Burgers equation = —2pt~^Vxx — 2pvVx studied in [17]. 

Note 3. The well-known Cole-Hopf transformation v = 2wx/w reduces equation (7) to the 
famous constant coefficient reaction-diffusion equation with weak nonlinearity 

wt = Wxx — 2pw In . 

After application of the Cole-Hopf transformation to the list of known exact solutions (see, 
e.g., [21]) of the equation with weak nonlinearity we obtain exactly solutions (8) of equation (7). 

In this and previous sections we have constructed exact solutions of variable coefficient equa- 
tions of form (1) reducible to constant coefficient ones which either belong or not to class (1). 
In the next section we consider an essentially variable coefficient equation that is distinguished 
with its symmetry properties from class (1). 

6 Group analysis of a remarkable s/ (2, ]R) -invariant 
diffusion— convection equation 

Analyzing the results of group classification, we can observe a number of G^-inequivalent equa- 
tions (1) which are invariant with respect to different realizations of the algebra ,s/(2,M). The 
set of such equations is practically exhausted by the well-known ( "constant coefficient" ) Burgers 
and n~^/^-diffusion equations and by the equations which are equivalent to them with respect 
to additional transformations (cases 3.16 and 3.15). This set is supplemented by the unique 
essentially variable coefficient equation 



(case 3.13). The si(2, M)-invariance of (9) is directly connected with the fact that h is not 
constant. The corresponding realization of the algebra sZ(2,M) is quite different from ones of 



and 4.3. 



X^Ut = {u ^^^Ux)x + x'^U. 



'X 



(9) 
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cases 3.16 and 3.15 and is the maximal Lie invariance algebra of equation (9). It was the reason 
for us to study equation (9) from the symmetry point of view in details. Indeed, instead of 
equation (9) we investigate the simpler equation 

x^vt = vv^x - -{vxf + x'^Vj; (10) 



for the function v = u~^^^, i.e. u = v~^^^. The maximal Lie invariance algebra yl™^^ of (10) is 
generated by the basis operators 

Pt = dt, D = tdt + xd^ + 3vd^, U = fdt + {2tx + x^)da, + 6{t + x)vdy. 

These operators generate the following one-parameter groups of point transformations: 

Pt : t = t + £, X = X, V = v; 

D: i = e^t, x = e^x, v = e^^v; 

Yi ~ ^ ~ t + X t _ V 



1-et' l-e{t + x) 1-et' {l-e{t + x))^' 

The complete Lie invariance group G^^^ is generated by both the above continuous transforma- 
tions and the discrete transformation of changing of sign in the triple {t, x, v). The transforma- 
tions from can be used for construction of new solutions from known ones. 

A list of proper G°^^-inequivalent subalgebras of is exhausted by the algebras {Pt), 

(D), (Pt + n), (Pt, D). 

The simplest Lie reduction of equation (10) to an algebraic equation is obtained with the 
algebra {Pt,D). Namely, the ansatz v = Cx^ results in the algebraic equation C(C — 2) = 
with respect to the constant C, i.e. either C = or C = 2. Thus, we have two solutions of 
equation (10) v = and v = 2x^ which are invariant with respect to the translation and scale 
transformations simultaneously. 

The above one-dimensional subalgebras reduce equation (10) to ordinary differential equa- 
tions. Let us list the corresponding ansatzes and reduced equations as well as some partial exact 
solutions of the reduced equations: 

(Pt) : V = ip{uj), uj = x, 
5 

(fPuju; - ^(fl = -Uj'^'Pu;; = C, (f = 2uj^ . (11) 

{D): v = t^ip{u), w = |, 
5 

f^Pojoj - ^(fil = 3w^</? - {uj + l)u}'^(pu,; (12) 
(^ = 0, ip = 2uj^, if = ^00^^ + 2oo^, ip = 2u;^{oj + lf, ip = ^oo'^iuo + if + 2(jj^{uj + if . 

(Pt -I- n) : V = {{t + x)"^ + l)^tp{ijj), Lo = arctan(f + x) — arctant, 
5 

ffPiOio - -^t = -fuj sin^w - 6(p'^. (13) 
b 

Here C is an arbitrary constant. 

Finally, we have the following set of G™^-inequivalent Lie invariant exact solutions of equa- 
tion (10) (below 5 G {0,1}): 

v = S, v = 2x^, v= — + 2x^, V = x{t + x)"^ l-^ + 2^ 
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We can extend the set of solutions with non-trivial invariance transformations: 

v = C{t + xf, v = 2^{t + xf, - = ^|^ + 2^(* + ^f> 

All the constructed solutions are polynomials with respect to x of degree not greater than 6 
with coefficients depending on t, i.e. they have the form 

6 

v = Y.^\t)x\ (14) 

i=0 

The set of all solutions of the form (14) is closed with respect to transformations from 

and is exhausted, up to translations with respect to t and scale transformations, by the above 

solutions w = (5, v = 5{t + x)^ and the solutions given by the generalized ansatz 

v = 2x^ + ^^{t)x^ + i/{t)x^ + /(^)a;^ (15) 

where ip^ solve the reduced system 

^t = 7^'--y)\ ^f = 18/-^¥.V, V^? = -^(V^^)^ + 2^V. (16) 

Reduction of equation (10) with ansatz (15) to system (16) is a consequence of generalized 
nonclassical invariance of (10) with respect to the third-order evolutionary operator 

{x^v^xx - 12x^Vj;a; + &^xv^ - 120u + I2x^)dy. 

See, e.g., [5,29] for a definition of generalized nonclassical (or conditional Lie-Backlund) sym- 
metries of evolution equations. 

System (16) can be reduced to the single third-order ordinary differential equation for the 
function f^: 

63ipfu + mv^f? + 126vjV4 + 192(9^^) V' + 16(v^')' = 0. (17) 

The set of all solutions of the form (15) is also closed with respect to transformations from 
and, therefore, naturally induces the invariance groups of system (16) and equation (17). 

Indeed, the maximal Lie invariance algebra A' of (17) is generated by the operators 

p' = dt, D' = tdt - ^"d^, , n' = t^dt + (6 - 2tp*)d^, . 

Using one-dimensional subalgebras of A' ^ we can reduce (17) to algebraic equations and, after 
solving them we construct Lie invariant solutions of (17) and, therefore, of (10) and (16). Thus, 
the subalgebra {P') leads to the simplest ansatz tp^ = C and the algebraic reduced equation 
C = 0. The subalgebra {P' + 11') gives the ansatz ip'^ = {6t + C)/{iP' + 1) and incompatible 
over M algebraic equation (4C^ + 81) (C^ -|- 36) = 0. The most interesting reduction of such 
type is connected with the subalgebra {D'). The (£)') -invariant ansatz (^^ = C/t reduces (17) 
to the algebraic equation 16^^ - 192C^ + 6390^ - 378C = having non-trivial set of solutions 
{0, 3/4, 21/4, 6}. The corresponding invariant solutions of equation (10) are adduced above. 

It is obvious that equation (9) has the nonclassical symmetry operator dt + dx which gives the 
ansatz v = (^(o;), u = t + x and reduces (9) to the equation (fif^juj — = 0- The corresponding 
solution V = Ci{x + t + Cq)^, where Co and Ci are arbitrary constants, is Lie invariant and 
equivalent to one adduced above. Let us note that the same operator dt -|- is a nonclassical 
symmetry of any equation (1), where f = h and g = 1. Detailed discussion of nonclassical 
(reduction) operators for equations (1) is given in the next section. 
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7 On nonclassical symmetries 



Reduction operators (nonclassical symmetries, Q-conditional symmetries) of equations (1) have 
the general form Q = rdf + ^d^ + r]du, where r, ^ and r] are functions of t, x and u, and 
(r,6 7^(0,0). 

Definition 2. The differential equation C is called conditionally invariant with respect to 
the operator Q if the relation a;, M(r))|^^g(^)= holds, which is called the conditional 

invariance criterion. Then Q is called an operator of conditional symmetry (or Q-conditional 
symmetry, nonclassical symmetry etc) of the equation jC. 

See, e.g., [26,31] for necessary definitions and properties of nonclassical symmetries. Since (1) 
is an evolution equation, there are two principally different cases of finding Q: t ^ and r = 0. 

The problem of classification of nonclassical symmetries with r 7^ is not completely solved 
even in the case of constant coefficient diffusion-convection equations. However, some particular 
results are known. Thus, e.g., constant coefficient diffusion equation with power nonlinearity 

admits non-classical symmetries only for /i = —1/2 [9]. Ibid a particular case of such operators 
Q = dt + 12x~'^v}/'^ is adduced and the corresponding exact solution u = (6ta;~^ + Cia,"^+C2X~^)^ 
is constructed. Here C\ and C2 are arbitrary constants. The same solution were found previously 
in [16]. Conditional symmetry operator Q = dt — Xu^^^dx for the equation ut = {u^Ux)x + 
\u^~^^Ux is found in the recent work [3]. Constant coefficient diffusion equation with exponential 
nonlinearity ut = {e^Ux)x admits two inequivalent reduction operators: Q = xdt + e"5u and 
Q = x'^dt + 2xe"5a; + 2e"5„ [6, 7]. 

We derive the system of determining equations for the reduction operators with r 7^ of 
equations (1): 

^uuA - ^uAu = (18) 
-VuuA^ - VuAAu + 2ixuA^ - 2iiufA - 2^uhAB - r^AA^u + r]{Auf = (19) 
-2r]xuA^ - 2r]xAAu - itfA + 2??e«/^ + ixxA^ - 2^^^/^ - ^^hAB 

ri^fAu + vhAuB - ff^A - ^h^AB - r^hAB^ = (20) 
r]tfA - A'^rixx - ri^hAB + 2^xVfA - rj^fA^ + ^rjf^A = (21) 

These equations can be partially integrated. In particular, (18) implies that 
^ = (t>{x,t)jA + i;{x,t). 

Then, (19) can be written in the form 



(^Mk^ =2{cPxA-cP^fJA-cP^Pf-cphB). 



In the case A = u'^ and B = (n 7^ — 1, — 2, — |, m 7^ — 1, — (n + 2)) the latter two formulas 
give 

e = 0(x,t)^."+i+V(x,O, 

+ -^Vn^+3_2(n + l)^^^^._ 



(n + 2)(2n + 3) ' n + 2 

Substituting these to (20) and (21) and solving the derived system we get an example of non- 
classical symmetry for equations (1). 
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Example 1. A = B = u V/(x), h{x), g{x) = 1: Q = dt + ^-y/uS-r, where (f) = (j){x) is an 
arbitrary solution of equation 

(j)" + h{x)<l)' -l/2f{x)(t)'^ = 

Similarly we found an example of nonclassical symmetry in the case of exponential nonlin- 
earities. 

Example 2. A = B = e^, dt + (j){x)e'^dx-i where ^' — fcjP' — hep = 0. 

As well-known, the operators with the vanishing coefficient of dt form so-called "no-go" case 
in study of conditional symmetries of an arbitrary (1 -|- l)-dimensional evolution equation since 
the problem on their finding is reduced to a single equation which is equivalent to the initial one 
(see e.g [8,22,30]). Note that "no-go" case has to be treated as impossibility only of exhaustive 
solving of the problem. A number of particular examples of reduction operators with r = can 
be constructed under additional constraints and then applied to finding exact solutions of the 
initial equation. Since the determining equation has more independent variables and, therefore, 
more degrees of freedom, it is more convenient often to guess a simple solution or a simple ansatz 
for the determining equation, which can give a parametric set of complicated solutions of the 
initial equation. Namely, in the case r = we have ^ / 0. Up to usual equivalence of reduction 
operators, ^ can be assumed equal to 1, i.e. Q = dx + rjdu- The conditional invariance criterion 
implies the determining equation on the coefficient r) 

7]tf - 27]7]uxfA - rj^'quufA - 2ifr]ufAu + rjrjufxA - rjxxfA - SrirjxfAu + rjxfxA 
^VxfhB - rffA^u + rffxAu + vfxhB - r^hxfB - iffhBu = 

which is reduced with a non-point transformation to equation (1), where r] becomes a parameter. 
We have found sonic partial solutions of the determining equations. 

Example 3. A = B, V/, V/i, r? = (t){x)A-^, where 4>' + hep = 0. 

_ 2(n+l) 2_ 

Example 4. A = u", B = u^^, f = x ^, h = x ■^+^ , rj = -r^W^ ^+^u'- 



8 Contractions of solutions 

Consider, as in Section 3, the class {C^} of totally nondegenerate systems C^: L{x,U(^p^, X) = 
of I differential equations for m unknown functions u = {u^, . . . ,u^) of n independent variables 
X = (xi, . . . which are parameterized with the (single numeric) parameter A. 

Let for any A from a deleted neighborhood of Aq the system C'^ be Lie invariant with respect 
to a symmetry algebra A^, ^ C^^ and A^ — > A^^ , A ^ Aq. Then, according to the results 
of Section 3, the system is Lie invariant with respect to the algebra A^° . 

If ui'^ is a Lie ansatz corresponding to the algebra then, using the theorems on convergence 
of solutions of differential equations with respect to initial data and parameters (see, e.g., [4]) one 
can prove that there exists such w^" that w'^ — > oo^° , A — >■ Aq and u)^^ is an Ansatz corresponding 
to the algebra A^° . Consequently the systems reduced with respect to converge to one reduced 
with Lo'^^ . In such sense we can talk about contractions of ansatzes and reduced systems. 

Consider an example of contractions of ansatzes, reduced equations and solutions of equa- 
tions (1). As it is shown above, under the action of contraction u = 1 -|- ^, /x — > -|-oo diffusion 
equation with power nonlinearity 

ut = {u''ux)x (22) 
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(case 3.14) invariant with respect to A^^ - 
xdx + 2^~^udu) is reduced to the equation 



(23) 



(case 2.7) with exponential nonUnearity being invariant with respect to the Lie symmetry algebra 
A^-p = (Qi = du Q2 = tdt - du, Qa = dx, Qa = xd^ + 2du). 

The basis elements of the Lie invariance algebra A^^^ can be obtained under the same con- 
traction as Q'' Qi. 

All possible inequivalent (with respect to inner automorphisms) one-dimensional subalgebras 
of the maximal Lie invariance algebras Af^ are exhausted by the ones listed in Table 5 together 
with the corresponding ansatzes and the reduced ODEs. 

Table 5. Reduced ODEs for (22). /x / 0, -1, a / 0, e = ±1, 5 = signt. 



N 



Subalgebra 



Ansatz u = 



t 
t 

X 
X 

X — et 
X — eln \t\ 
xe-^' 
x\t\-'' 



Reduced ODE 



(p' = 

V3' = 2/i-2(2 + /i)^'^+i 



= 

= W 



= —euji{>' + 2/Lt ^£ifi 



{ifi^ip'y = 5ii ^ {2a — V)ifi — Sauap 



For the considered equations the optimal systems of subalgebras of A^^ are contracted to an 
optimal system of subalgebras of A^^"^. Let us note, that in general the question whether optimal 
systems of subalgebras of the maximal Lie invariance algebras of systems converge to the 
optimal subalgebras system of >C^°, remains open. 

Contractions of ansatzes and reduced ODEs in cases 5.1 and 5.3 are obvious. Let us consider 
in more details contraction of case 5.4. Under the transformation of variables u = 1 + the 



function should be changed as 
follows: 



1 -)- ^. The ansatz u 
1^ 



ip{x)\t\ can be contracted as 



(i+5r 



i + ^)V 



s" = em 



\i — > oo, 



Therefore, u = p> — \x\\t\. The similarity variable lj = x is not changed under the contraction. 
Substituting the derived expressions to the reduced equation 4.4 we obtain 



M + 1 



-(5 1 + 



If now /i ^ oo we get the reduced ordinary differential equation 
{f)" = -S 

for the target equation (23). 

Similarly one can contract all ansatzes and reduced equations of (22) to ones for equations 
with exponential nonlinear ity. The results of these contractions are summarized in Table 6. 
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Table 6. Reduced ODEs for (23). a^O, s = ±1, S = signt. 



N 


Subalgebra 


Ansatz u = 


to 


Reduced ODE 


1 






t 


= 


2 




(p{uj) + 21n \x\ 


t 


ip' = 26"^ 


3 


(Qi) 


ip{ijj) 


X 


(e^)" = 


4 


(Q2> 


ffiiio) — In \t\ 


X 


(e^)" = -S 


5 


(Qi+eQa) 




X — et 


{en" = -e^' 


6 


{Q2 + eQa) 


fiuj) — In \t\ 


a; — eln \t\ 


(ef)" = -5{e^' + 1) 


7 


{Qi+eQi) 


(p{uj) + 2et 


xe-^^ 


(e^)" = -ecjifi' + 2e 


8 


{Q2 + aQi) 


ifiui) + (2a - 1) In \t\ 


xW^ 


{e^y = S{-aLOifi' + 2a - 1) 



9 Conclusion 

In this second part of the presented series of papers (see also [12-14]) we investigate in more 
detail symmetry properties of class (1). Namely, considering non-trivial limits of parameterized 
subclasses of equations from class (1), which generate contractions of the corresponding maximal 
Lie invariance algebras, we introduce the notion of contraction of (systems of) differential equa- 
tions and consider contractions of equations and from class (1) and ones of their symmetries. 
We also investigate si(2, M)-invariant equation (9) which is "essentially variable coefficient" in 
the sense that it is not reducible to equations of form (1) with constant values of /, g and h. 

Using similar techniques, we can study other classes of non-linear evolution equations which 
arc closed to the class under consideration, e.g. the class of variable coefficient reaction-diffusion 
equations of the general form f{x)ut = {g{x)A{u)ux)x+h{x)B(u), where all denotations coincide 
with ones in (1). (See the resent work [27] for its particular case with power nonlinearity.) 
However, experience of modern group analysis shows that extension of circle of problems leads 
to necessity of modification and enhancement of applied tools. 

Another natural direction is generalization of the "contraction concept" to conservation laws 
of (systems of) diffusion equations that is considered in the next part [13] of this series. 
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